Personalized neoantigen vaccine and pembrolizumab in advanced hepatocellular carcinoma: a phase 1/2 trial.

Programmed cell death protein 1 (PD-1) inhibitors have modest efficacy as a monotherapy in hepatocellular carcinoma (HCC). A personalized therapeutic cancer vaccine (PTCV) may enhance responses to PD-1 inhibitors through the induction of tumor-specific immunity. We present results from a single-arm, open-label, phase 1/2 study of a DNA plasmid PTCV (GNOS-PV02) encoding up to 40 neoantigens coadministered with plasmid-encoded interleukin-12 plus pembrolizumab in patients with advanced HCC previously treated with a multityrosine kinase inhibitor. Safety and immunogenicity were assessed as primary endpoints, and treatment efficacy and feasibility were evaluated as secondary endpoints. The most common treatment-related adverse events were injection-site reactions, observed in 15 of 36 (41.6%) patients. No dose-limiting toxicities or treatment-related grade ≥3 events were observed. The objective response rate (modified intention-to-treat) per Response Evaluation Criteria in Solid Tumors 1.1 was 30.6% (11 of 36 patients), with 8.3% (3 of 36) of patients achieving a complete response. Clinical responses were associated with the number of neoantigens encoded in the vaccine. Neoantigen-specific T cell responses were confirmed in 19 of 22 (86.4%) evaluable patients by enzyme-linked immunosorbent spot assays. Multiparametric cellular profiling revealed active, proliferative and cytolytic vaccine-specific CD4+ and CD8+ effector T cells. T cell receptor ß-chain (TCRß) bulk sequencing results demonstrated vaccination-enriched T cell clone expansion and tumor infiltration. Single-cell analysis revealed posttreatment T cell clonal expansion of cytotoxic T cell phenotypes. TCR complementarity-determining region cloning of expanded T cell clones in the tumors following vaccination confirmed reactivity against vaccine-encoded neoantigens. Our results support the PTCV's mechanism of action based on the induction of antitumor T cells and show that a PTCV plus pembrolizumab has clinical activity in advanced HCC. ClinicalTrials.gov identifier: NCT04251117 .

Analytical validation of NeXT Personal®, an ultra-sensitive personalized circulating tumor DNA assay.

We describe the analytical validation of NeXT Personal®, an ultra-sensitive, tumor-informed circulating tumor DNA (ctDNA) assay for detecting residual disease, monitoring therapy response, and detecting recurrence in patients diagnosed with solid tumor cancers. NeXT Personal uses whole genome sequencing of tumor and matched normal samples combined with advanced analytics to accurately identify up to ~1,800 somatic variants specific to the patient's tumor. A personalized panel is created, targeting these variants and then used to sequence cell-free DNA extracted from patient plasma samples for ultra-sensitive detection of ctDNA. The NeXT Personal analytical validation is based on panels designed from tumor and matched normal samples from two cell lines, and from 123 patients across nine cancer types. Analytical measurements demonstrated a detection threshold of 1.67 parts per million (PPM) with a limit of detection at 95% (LOD95) of 3.45 PPM. NeXT Personal showed linearity over a range of 0.8 to 300,000 PPM (Pearson correlation coefficient = 0.9998). Precision varied from a coefficient of variation of 12.8% to 3.6% over a range of 25 to 25,000 PPM. The assay targets 99.9% specificity, with this validation study measuring 100% specificity and in silico methods giving us a confidence interval of 99.92 to 100%. In summary, this study demonstrates NeXT Personal as an ultra-sensitive, highly quantitative and robust ctDNA assay that can be used to detect residual disease, monitor treatment response, and detect recurrence in patients.

Analytic validation of NeXT Dx™, a comprehensive genomic profiling assay.

We describe the analytic validation of NeXT Dx, a comprehensive genomic profiling assay to aid therapy and clinical trial selection for patients diagnosed with solid tumor cancers. Proprietary methods were utilized to perform whole exome and whole transcriptome sequencing for detection of single nucleotide variants (SNVs), insertions/deletions (indels), copy number alterations (CNAs), and gene fusions, and determination of tumor mutation burden and microsatellite instability. Variant calling is enhanced by sequencing a patient-specific normal sample from, for example, a blood specimen. This provides highly accurate somatic variant calls as well as the incidental reporting of pathogenic and likely pathogenic germline alterations. Fusion detection via RNA sequencing provides more extensive and accurate fusion calling compared to DNA-based tests. NeXT Dx features the proprietary Accuracy and Content Enhanced technology, developed to optimize sequencing and provide more uniform coverage across the exome. The exome was validated at a median sequencing depth of >500x. While variants from 401 cancer-associated genes are currently reported from the assay, the exome/transcriptome assay is broadly validated to enable reporting of additional variants as they become clinically relevant. NeXT Dx demonstrated analytic sensitivities as follows: SNVs (99.4%), indels (98.2%), CNAs (98.0%), and fusions (95.8%). The overall analytic specificity was >99.0%.

Precision Neoantigen Discovery Using Large-Scale Immunopeptidomes and Composite Modeling of MHC Peptide Presentation.

Major histocompatibility complex (MHC)-bound peptides that originate from tumor-specific genetic alterations, known as neoantigens, are an important class of anticancer therapeutic targets. Accurately predicting peptide presentation by MHC complexes is a key aspect of discovering therapeutically relevant neoantigens. Technological improvements in mass spectrometry-based immunopeptidomics and advanced modeling techniques have vastly improved MHC presentation prediction over the past 2 decades. However, improvement in the accuracy of prediction algorithms is needed for clinical applications like the development of personalized cancer vaccines, the discovery of biomarkers for response to immunotherapies, and the quantification of autoimmune risk in gene therapies. Toward this end, we generated allele-specific immunopeptidomics data using 25 monoallelic cell lines and created Systematic Human Leukocyte Antigen (HLA) Epitope Ranking Pan Algorithm (SHERPA), a pan-allelic MHC-peptide algorithm for predicting MHC-peptide binding and presentation. In contrast to previously published large-scale monoallelic data, we used an HLA-null K562 parental cell line and a stable transfection of HLA allele to better emulate native presentation. Our dataset includes five previously unprofiled alleles that expand MHC diversity in the training data and extend allelic coverage in underprofiled populations. To improve generalizability, SHERPA systematically integrates 128 monoallelic and 384 multiallelic samples with publicly available immunoproteomics data and binding assay data. Using this dataset, we developed two features that empirically estimate the propensities of genes and specific regions within gene bodies to engender immunopeptides to represent antigen processing. Using a composite model constructed with gradient boosting decision trees, multiallelic deconvolution, and 2.15 million peptides encompassing 167 alleles, we achieved a 1.44-fold improvement of positive predictive value compared with existing tools when evaluated on independent monoallelic datasets and a 1.17-fold improvement when evaluating on tumor samples. With a high degree of accuracy, SHERPA has the potential to enable precision neoantigen discovery for future clinical applications.

A machine learning algorithm with subclonal sensitivity reveals widespread pan-cancer human leukocyte antigen loss of heterozygosity.

Human leukocyte antigen loss of heterozygosity (HLA LOH) allows cancer cells to escape immune recognition by deleting HLA alleles, causing the suppressed presentation of tumor neoantigens. Despite its importance in immunotherapy response, few methods exist to detect HLA LOH, and their accuracy is not well understood. Here, we develop DASH (Deletion of Allele-Specific HLAs), a machine learning-based algorithm to detect HLA LOH from paired tumor-normal sequencing data. With cell line mixtures, we demonstrate increased sensitivity compared to previously published tools. Moreover, our patient-specific digital PCR validation approach provides a sensitive, robust orthogonal approach that could be used for clinical validation. Using DASH on 610 patients across 15 tumor types, we find that 18% of patients have HLA LOH. Moreover, we show inflated HLA LOH rates compared to genome-wide LOH and correlations between CD274 (encodes PD-L1) expression and microsatellite instability status, suggesting the HLA LOH is a key immune resistance strategy.

Prediction of Immunotherapy Response in Melanoma through Combined Modeling of Neoantigen Burden and Immune-Related Resistance Mechanisms.

PURPOSE: While immune checkpoint blockade (ICB) has become a pillar of cancer treatment, biomarkers that consistently predict patient response remain elusive due to the complex mechanisms driving immune response to tumors. We hypothesized that a multi-dimensional approach modeling both tumor and immune-related molecular mechanisms would better predict ICB response than simpler mutation-focused biomarkers, such as tumor mutational burden (TMB). EXPERIMENTAL DESIGN: Tumors from a cohort of patients with late-stage melanoma (n = 51) were profiled using an immune-enhanced exome and transcriptome platform. We demonstrate increasing predictive power with deeper modeling of neoantigens and immune-related resistance mechanisms to ICB. RESULTS: Our neoantigen burden score, which integrates both exome and transcriptome features, more significantly stratified responders and nonresponders (P = 0.016) than TMB alone (P = 0.049). Extension of this model to include immune-related resistance mechanisms affecting the antigen presentation machinery, such as HLA allele-specific LOH, resulted in a composite neoantigen presentation score (NEOPS) that demonstrated further increased association with therapy response (P = 0.002). CONCLUSIONS: NEOPS proved the statistically strongest biomarker compared with all single-gene biomarkers, expression signatures, and TMB biomarkers evaluated in this cohort. Subsequent confirmation of these findings in an independent cohort of patients (n = 110) suggests that NEOPS is a robust, novel biomarker of ICB response in melanoma.

Precision Neoantigen Discovery Using Large-scale Immunopeptidomes and Composite Modeling of MHC Peptide Presentation.

Major histocompatibility complex (MHC)-bound peptides that originate from tumor-specific genetic alterations, known as neoantigens, are an important class of anticancer therapeutic targets. Accurately predicting peptide presentation by MHC complexes is a key aspect of discovering therapeutically relevant neoantigens. Technological improvements in mass-spectrometry-based immunopeptidomics and advanced modeling techniques have vastly improved MHC presentation prediction over the past two decades. However, improvement in the sensitivity and specificity of prediction algorithms is needed for clinical applications such as the development of personalized cancer vaccines, the discovery of biomarkers for response to checkpoint blockade, and the quantification of autoimmune risk in gene therapies. Toward this end, we generated allele-specific immunopeptidomics data using 25 monoallelic cell lines and created Systematic HLA Epitope Ranking Pan Algorithm (SHERPA), a pan-allelic MHC-peptide algorithm for predicting MHC-peptide binding and presentation. In contrast to previously published large-scale monoallelic data, we used an HLA-null K562 parental cell line and a stable transfection of HLA alleles to better emulate native presentation. Our dataset includes five previously unprofiled alleles that expand MHC-binding pocket diversity in the training data and extend allelic coverage in under profiled populations. To improve generalizability, SHERPA systematically integrates 128 monoallelic and 384 multiallelic samples with publicly available immunoproteomics data and binding assay data. Using this dataset, we developed two features that empirically estimate the propensities of genes and specific regions within gene bodies to engender immunopeptides to represent antigen processing. Using a composite model constructed with gradient boosting decision trees, multiallelic deconvolution, and 2.15 million peptides encompassing 167 alleles, we achieved a 1.44-fold improvement of positive predictive value compared with existing tools when evaluated on independent monoallelic datasets and a 1.15-fold improvement when evaluating on tumor samples. With a high degree of accuracy, SHERPA has the potential to enable precision neoantigen discovery for future clinical applications.

High prevalence of focal and multi-focal somatic genetic variants in the human brain.

Somatic mutations during stem cell division are responsible for several cancers. In principle, a similar process could occur during the intense cell proliferation accompanying human brain development, leading to the accumulation of regionally distributed foci of mutations. Using dual platform >5000-fold depth sequencing of 102 genes in 173 adult human brain samples, we detect and validate somatic mutations in 27 of 54 brains. Using a mathematical model of neurodevelopment and approximate Bayesian inference, we predict that macroscopic islands of pathologically mutated neurons are likely to be common in the general population. The detected mutation spectrum also includes DNMT3A and TET2 which are likely to have originated from blood cell lineages. Together, these findings establish developmental mutagenesis as a potential mechanism for neurodegenerative disorders, and provide a novel mechanism for the regional onset and focal pathology in sporadic cases.

svclassify: a method to establish benchmark structural variant calls.

BACKGROUND: The human genome contains variants ranging in size from small single nucleotide polymorphisms (SNPs) to large structural variants (SVs). High-quality benchmark small variant calls for the pilot National Institute of Standards and Technology (NIST) Reference Material (NA12878) have been developed by the Genome in a Bottle Consortium, but no similar high-quality benchmark SV calls exist for this genome. Since SV callers output highly discordant results, we developed methods to combine multiple forms of evidence from multiple sequencing technologies to classify candidate SVs into likely true or false positives. Our method (svclassify) calculates annotations from one or more aligned bam files from many high-throughput sequencing technologies, and then builds a one-class model using these annotations to classify candidate SVs as likely true or false positives. RESULTS: We first used pedigree analysis to develop a set of high-confidence breakpoint-resolved large deletions. We then used svclassify to cluster and classify these deletions as well as a set of high-confidence deletions from the 1000 Genomes Project and a set of breakpoint-resolved complex insertions from Spiral Genetics. We find that likely SVs cluster separately from likely non-SVs based on our annotations, and that the SVs cluster into different types of deletions. We then developed a supervised one-class classification method that uses a training set of random non-SV regions to determine whether candidate SVs have abnormal annotations different from most of the genome. To test this classification method, we use our pedigree-based breakpoint-resolved SVs, SVs validated by the 1000 Genomes Project, and assembly-based breakpoint-resolved insertions, along with semi-automated visualization using svviz. CONCLUSIONS: We find that candidate SVs with high scores from multiple technologies have high concordance with PCR validation and an orthogonal consensus method MetaSV (99.7 % concordant), and candidate SVs with low scores are questionable. We distribute a set of 2676 high-confidence deletions and 68 high-confidence insertions with high svclassify scores from these call sets for benchmarking SV callers. We expect these methods to be particularly useful for establishing high-confidence SV calls for benchmark samples that have been characterized by multiple technologies.

Achieving high-sensitivity for clinical applications using augmented exome sequencing.

BACKGROUND: Whole exome sequencing is increasingly used for the clinical evaluation of genetic disease, yet the variation of coverage and sensitivity over medically relevant parts of the genome remains poorly understood. Several sequencing-based assays continue to provide coverage that is inadequate for clinical assessment. METHODS: Using sequence data obtained from the NA12878 reference sample and pre-defined lists of medically-relevant protein-coding and noncoding sequences, we compared the breadth and depth of coverage obtained among four commercial exome capture platforms and whole genome sequencing. In addition, we evaluated the performance of an augmented exome strategy, ACE, that extends coverage in medically relevant regions and enhances coverage in areas that are challenging to sequence. Leveraging reference call-sets, we also examined the effects of improved coverage on variant detection sensitivity. RESULTS: We observed coverage shortfalls with each of the conventional exome-capture and whole-genome platforms across several medically interpretable genes. These gaps included areas of the genome required for reporting recently established secondary findings (ACMG) and known disease-associated loci. The augmented exome strategy recovered many of these gaps, resulting in improved coverage in these areas. At clinically-relevant coverage levels (100 % bases covered at ≥20×), ACE improved coverage among genes in the medically interpretable genome (>90 % covered relative to 10-78 % with other platforms), the set of ACMG secondary finding genes (91 % covered relative to 4-75 % with other platforms) and a subset of variants known to be associated with human disease (99 % covered relative to 52-95 % with other platforms). Improved coverage translated into improvements in sensitivity, with ACE variant detection sensitivities (>97.5 % SNVs, >92.5 % InDels) exceeding that observed with conventional whole-exome and whole-genome platforms. CONCLUSIONS: Clinicians should consider analytical performance when making clinical assessments, given that even a few missed variants can lead to reporting false negative results. An augmented exome strategy provides a level of coverage not achievable with other platforms, thus addressing concerns regarding the lack of sensitivity in clinically important regions. In clinical applications where comprehensive coverage of medically interpretable areas of the genome requires higher localized sequencing depth, an augmented exome approach offers both cost and performance advantages over other sequencing-based tests.

Variant priorization and analysis incorporating problematic regions of the genome.

In case-control studies of rare Mendelian disorders and complex diseases, the power to detect variant and gene-level associations of a given effect size is limited by the size of the study sample. Paradoxically, low statistical power may increase the likelihood that a statistically significant finding is also a false positive. The prioritization of variants based on call quality, putative effects on protein function, the predicted degree of deleteriousness, and allele frequency is often used as a mechanism for reducing the occurrence of false positives, while preserving the set of variants most likely to contain true disease associations. We propose that specificity can be further improved by considering errors that are specific to the regions of the genome being sequenced. These problematic regions (PRs) are identified a-priori and are used to down-weight constitutive variants in a case-control analysis. Using samples drawn from 1000-Genomes, we illustrate the utility of PRs in identifying true variant and gene associations using a case-control study on a known Mendelian disease, cystic fibrosis (CF).

Methane biogenesis during sodium azide-induced chemical hypoxia in rats.

Previous studies demonstrated methane generation in aerobic cells. Our aims were to investigate the methanogenic features of sodium azide (NaN(3))-induced chemical hypoxia in the whole animal and to study the effects of l-α-glycerylphosphorylcholine (GPC) on endogenous methane production and inflammatory events as indicators of a NaN(3)-elicited mitochondrial dysfunction. Group 1 of Sprague-Dawley rats served as the sham-operated control; in group 2, the animals were treated with NaN(3) (14 mg·kg(-1)·day(-1) sc) for 8 days. In group 3, the chronic NaN(3) administration was supplemented with daily oral GPC treatment. Group 4 served as an oral antibiotic-treated control (rifaximin, 10 mg·kg(-1)·day(-1)) targeting the intestinal bacterial flora, while group 5 received this antibiotic in parallel with NaN(3) treatment. The whole body methane production of the rats was measured by means of a newly developed method based on photoacoustic spectroscopy, the microcirculation of the liver was observed by intravital videomicroscopy, and structural changes were assessed via in vivo fluorescent confocal laser-scanning microscopy. NaN(3) administration induced a significant inflammatory reaction and methane generation independently of the methanogenic flora. After 8 days, the hepatic microcirculation was disturbed and the ATP content was decreased, without major structural damage. Methane generation, the hepatic microcirculatory changes, and the increased tissue myeloperoxidase and xanthine oxidoreductase activities were reduced by GPC treatment. In conclusion, the results suggest that methane production in mammals is connected with hypoxic events associated with a mitochondrial dysfunction. GPC is protective against the inflammatory consequences of a hypoxic reaction that might involve cellular or mitochondrial methane generation.

Retroviral activation of the mir-106a microRNA cistron in T lymphoma.

Retroviral insertion into a host genome is a powerful tool not only for the discovery of cancer genes, but also for the discovery of potential oncogenic noncoding RNAs. In a large-scale mouse T lymphocyte tumor screen we found a high density of integrations upstream of the mir-106a microRNA cistron. In tumors containing an integration, the primary transcript encoding the mir-106a cistron was overexpressed five to 20-fold compared with that of control tumors; concomitantly, the mature mir-106a and mir-363 microRNAs were highly overexpressed as well. These findings suggest the mir-106a cistron plays an important role in T cell tumorigenesis.

Activation of an oncogenic microRNA cistron by provirus integration.

Retroviruses can cause tumors when they integrate near a protooncogene or tumor suppressor gene of the host. We infected >2,500 mice with the SL3-3 murine leukemia virus; in 22 resulting tumors, we found provirus integrations nearby or within the gene that contains the mir-17-92 microRNA (miRNA) cistron. Using quantitative real-time PCR, we showed that expression of miRNA was increased in these tumors, indicating that retroviral infection can induce expression of oncogenic miRNAs. Our results demonstrate that retroviral mutagenesis can be a potent tool for miRNA discovery.